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The results of pion fluctuation measurements in SERP-E-190 experiment (project Ther- 
malization) with 50 GeV proton beam irradiation of the liquid hydrogen target at SVD-2 
setup are presented. The photons are detected in the electromagnetic calorimeter. MC mod- 
cling of photon detection has shown the linear dependence between number of photons in the 
calorimeter and the average number of neutral pions. Neutral pion number Nq distributions 
for each total number of particles in an event N to t = N c h + Nq are obtained after making 
corrections on the setup acceptance, triggering and efficiency of the event reconstruction. 
The scaled variance of neutral pion fluctuations, ui = D / < No >, is measured. The fluc- 
tuations increase at N to t >22. According to quantum statistics models it may indicate for 
the approaching to pion condensate conditions for high pion multiplicity in the system. This 
effect have been observed for the first time. 
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I. INTRODUCTION 

The SERP-E-190 experiment (project Ther- 
malization) is carried out at upgraded setup 
SVD-2 (Spectrometer with Vertex Detector), 
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which is irradiated with 50 GeV proton beam 
from U-70 IHEP (Protvino) accelerator. The 
setup is equipped with a liquid hydrogen target, 
a microstrip silicon vertex detector (VD), a drift 
tube tracker, magnetic spectrometer, Cherenkov 
counter and an electromagnetic calorimeter for 
photon detection (DEGA) [1]. Project Thermal- 
ization is aimed at studying processes with high 
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multiplicity in pp- interactions, which are one of 
fundamental regions of hadron physics. It is 
not possible to describe them with perturbative 
QCD. The theory gives only a qualitative pic- 
ture of these processes. The multiplicity distri- 
bution at 50 GeV was measured earlier [2] up to 
the number of charged particles N c h=16 with av- 
erage number of charged particles at this energy 
< N c h >=5.3. The kinematical limit for the total 
number charged and neutral particles is N tot =56. 
The data for N ch = 4^-22 and for N tot = 4-5-31 
are presented in this work. Collective effects can 
occur in the events with the multiplicity by sev- 
eral times higher than the average one: the large 
fluctuations of charged and neutral pion numbers, 
as result of formation pion condensate, formation 
of jets with identical pions, the so-called Bose- 
Einstein multiparticle effect, formation of ring 
events as result gluon hadronization, which are 
emitted out by partons in the nuclear environ- 
ment (Cherenkov radiation analog) and others. 
The SVD-2 data allows one to check and develop 
various models of multiparticle production with 

N tot >< N tot >■ 

M. I. Gorenstein and V. V. Begun [3] [4] have 
shown that at the approach of the pion system to 
Bose-Einstein condensate conditions (BEC) the 
neutral pion number fluctuations are increasing 
in accordance with the model based on quantum 
statistics. These fluctuations can be detected by 
the scaled variance, oo, which is defined as the 
ratio of variance D for neutral pion number dis- 
tribution Nq and average < iVo >> 

u = D/ < N >. 

The value of u rising with increasing of the total 
particle number, N to t = N c h + No, depends on 
temperature and energy density of pion system. 

A part of E-190 data is used in this paper 
to search for pion number fluctuations. Numbers 
of charged particles, N c h, and photons, iV 7 , in 
each event are detected. These values are cor- 
rected on the setup acceptance and reconstruc- 
tion efficiency by means of modeling. Simula- 
tion has allowed getting the neutral pion num- 
ber, No, in each event also. Numbers of events, 
N ev (No, N to t), are measured here with corrections 
on various losses. For the analysis of the data at 
different Ntot relative values of no = No/Ntot and 
n) = N ev (N , N tot ) /N ev (Ntot) are used. Thus n 
is changed in the range of 0-5-1 and the sum of 




FIG. 1: Distributions r$ for normalized multiplicity 
of neutral pions in QCD model and when the system 
approaches to BEC. 

all tq is equal to 1 for each N to t (normalization 
condition). The Fig. 1 qualitatively illustrates 
the behaviour of value ro for three cases: a) the 
generation of events with PYTHIA5.6 program, 
b) the pion system in which some pions drop out 
into condensate, c) all pions are in the BEC con- 
dition. Each distribution is characterized by av- 
erage, < no >, and by standard deviation, a for 
a Gaussian fit. 



II. SIMULATION OF NEUTRAL PION 
DETECTION 

Calorimeter DEGA at SVD-2 setup detects 
the events with photons from neutral pions de- 
cay. Registration of all 7r° in the event is not pos- 
sible because of limited DEGA aperture and the 
threshold on photon detection energy. But 7r° re- 
construction efficiency can be estimated by means 
of simulation. Using PYTHIA5.6 code 10 6 events 
(MC) are simulated for pp— >-X inelastic interac- 
tions at 50 GeV. Only the events with N c h >4 are 
analyzed. The photon detection efficiency is as- 
sumed to be equal to 1 if photon hits DEGA and 
its energy is greater than 100 MeV. Practically 
all photons are the product of 7r° decays (95%), 
37% of them give two gammas signal in DEGA 
and 18% of vr° result in one photon signal. Fig. 
2 shows the dependence of pion number in event, 
No, on number of photons in DEGA, N T It is 
clear that there is no unique connection between 
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FIG. 2: a) The dependence of the pions number in 
MC events on photons number detected in DEGA, 
b) < No > and N 7 correlation 



iV 7 and No. (Fig. 2a). Instead each iV 7 is as- 
sociated with some number of A"o and there is a 
linear correlation between average < Nq > and 
A" 7 (Fig. 2b). So relation between the number 
of events N ev (N^, N c h) and N ev (No, N c h) can be 
found from this analysis. In Fig. 2c multiplic- 
ity distributions for pions and photons in DEGA 
are presented. The form of these distributions is 
similar except the area of small multiplicity. 

Values n and r are calculated for MC events. 
Then the dependence of r versus n has been fitted 
with Gaussian. In Fig. 3a the fitted parameters 
< tiq > and < n 7 > on N to t are presented, where 
Ntot = N ch + Nq for 7T° and N tot = N ch + Nj for 
photons. Standard deviations, a of Gaussians, 
are presented in Fig. 3b too. The normalized 
dispersion, u=a 2 * N to t/ < n > is shown in Fig. 
3c. The value oj decreases for photons, but re- 
mains near the constant for pions in the full area 
of Ntot changes. 



III. PHOTON AND CHARGED 
PARTICLES RECONSTRUCTION 

The calorimeter DEGA contains 
42x32=1344 elements from lead glass blocks 
with PM. The calibration of DEGA is carried 
out with irradiation of each element center by 
15 GeV electron beam. Almost all (98%) energy 
of electromagnetic shower from the photon is 
located in the cell (3x3) elements. Photon 
reconstruction consists in the searching of (3x3) 
signal clusters and analyzing them with criteria 
for photon. The present work comprises ~ 5 * 10 5 
events of pp interaction. The analysis of the 
electromagnetic showers has led to the following 
results: average photon number in the event 
is < A^ 7 >=3.0 (Fig. 4), their average energy 
is < E 1 >=2.8 GeV (Fig. 5b), the minimum 
energy of photon detection is equal to 100 MeV. 

For the charged tracks reconstruction the 
data only from VD have been used. Because of 
various losses the corrections of charged particle 
number are essential for measurements of pion 
fluctuations. The correction for the setup accep- 
tance and the particle reconstruction efficiency 
is made by means of the apparatus performance 
simulation. The obtained weights have shown the 
contribution of events with a different true num- 
ber of charged tracks into the event sample with 
the reconstructed multiplicity. The details of this 
procedure are presented in [5] .In the present work 
the number of events with measured N c h is dis- 
tributed on event samples with restored value N c ^ 
according to these weights. 

Topological cross sections for pp- 
interactions at 50 GeV were measured in 
experiment with Mirabelle bubble chamber [2]. 
Taking into account that data topological cross 
sections for 10 < N c h < 24 had been obtained 
in [5]. The present data are obtained with 
suppression of small charged multiplicity events 
(N c h < 8) with trigger conditions. Therefore 
the charged multiplicity distribution has been 
corrected to lead to the distribution received in 
[2] [5]. 

Fig. 5 shows charged multiplicity distribu- 
tions after corrections. It is necessary to stress 
that the change of the event number for A^ af- 
ter introduction of corrections also leads to the 
change of the event number for N^. So, the dis- 
tribution for A~ 7 in Fig. 5 differs from the similar 
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FIG. 4: a) Photon multiplicity distribution; b) photon energy distribution for MC events. 
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FIG. 5: Distributions of corrected N c h, N 7 and N tot = 
N ch + N 7 . 

distribution in Fig. 4a for observable event num- 
bers. 



IV. NEUTRAL PION FLUCTUATION 
MEASUREMENTS 

Thus we have corrected event num- 
bers N ev (N J ,N c f l ). Then two-dimensional 
N ev (Nj, Nq) distributions for MC events (see 
Fig. 2) are used to recover event numbers 
N ev (No,N c h). We have introduced notations 
i=N 7 , i=N and N^N^Nq) = N ev (i,j). 
For each N c y t matrix of coefficients Cij = 
N e v(i,j)/N ev (i) is calculated, where N ev (i) = 
Y,j N ev{i,j)- Event number N ev (N~ p N ch ) is de- 
composed in sums of events with various Nq, 
N ev (i,j) = * N ev (i) at A r c / l =const. Normal- 
ization condition X^ c «j = l * s satisfied. Resulting 
sum N ev (j) = J2i N ev(i,j) at N ch =hx is the ana- 
log of event number N ev (N^, N c h), but for pions. 

The simulation by PYTHIA5.6 allows to ob- 
tain for V 7 < 10 and N c h < 14 only because 
of limitation of the MC events statistics. Regu- 
larities of factors C{j &X6 used to continue them 
to N-y > 10 and N c h > 14 region. Fig. 6a illus- 
trates the dependence of c^ factors on No for var- 
ious A" 7 and N^. The form of these distributions 
slightly depends on iV 7 and N c h, but their average 
< No > increases with 7V 7 . The dependence of 
the average < Nq > and standard deviation (rms) 



on iV 7 are shown in Fig. 6b. After fitting it by 
linear dependence the coefficients c^ for A" 7 >10 
and N c h >14 are calculated and the full sample 
of N ev (N to t, N c h, No) is obtained, which is used 
then to determine pion fluctuation. 

As mention before we have used scaled vari- 
ables no and ro (see Introduction): no = No/N tot 
and r (n ) = N ev (N , N tot )/N ev (N tot ), where 
Ntot = Nq + iV c /j. Function ro(no) are shown in 
Fig. 7 for every N tot > 10. All distributions are 
fitted with Gaussian function. The dependence 
of the fitting parameters is presented in Fig. 8. 
The data in the intervals (26, 27, 28) and (29, 30, 
31) are combined due to small statistics. 

One can see that the measured average < 
no > (Fig. 8) coincide with the same values for 
the neutral pions from MC events at Ntot >12. In 
the gluon dominance model (MGD) [6] neutral pi- 
ons average dependence on N to t has been received 
by analytical way. This dependence is also pre- 
sented in Fig. 8 and illustrates quite good agree- 
ment with the experimental data at N to t >14. 
The average < n 7 > is also shown. The measured 
standard deviations, a (Gaussian), (Fig. 8b) have 
shown the qualitative agreement with MC model 
only for N to t <22. The measured values a in- 
crease at higher N tot . 

The theoretical prediction of scaled variance 
ijj behavior (in our case uj = D(No)/ < No >= 
a 2 * Ntot/ < no >) is given in [3]. The analysis 
has been done for three energy densities of the 
pion system at the approach to the Bose-Einstein 
condensate condition (pion condensate) (Fig. 9). 
Our experimental data (Fig. 9b) have confirmed 
assumption on the BEC formation in pion system 
at Ntot >22 in pp-interactions at 50 GeV. 



V. CONCLUSION 

Measurements of the charged and neutral pi- 
ons number in the events with high multiplicity 
in pp-interactions at 50 GeV (experiment SERP— 
190, SVD-2 setup) together with MC analysis led 
to the following results: 

• The average number of neutral pions in the 
event is proportional to the photons num- 
ber detected in DEGA calorimeter that al- 
lows one to extract pion number fluctua- 
tions from photon number fluctuations. 



5 



SVD-2 Collaboration 




FIG. 6: a) The dependence of dj factors on iVo for various N 1 and N c h, b) fitting parameters: < No > and 
standard deviation (rms) as function of iV 7 for various N c h- 




FIG. 7: Scaled neutral pions 
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• It is convenient to present data in the 
scaled form: uq = N^/Ntot and ^o(no) = 
N ev (N , Ntot) /N ev {Ntot) with interval n is 
equal to CM-1. 

• The corrections for limited aperture VD, 
trigger action and efficiency of data pro- 
cessing system have been introduced to the 
data. 

• The function ro(rio) is fitted with Gaus- 
sian function and values < uq >, a and 
scaled variance u = D/ < Nq > are de- 
rived. They have shown the qualitative 



agreement with the same values received 
for PYTHIA5.6 code at N tot <22. 

• Pion number fluctuations increase at 
Ntot >22, that indicates approaching to 
pion condensate conditions for the high 
multiplicity pion system according to GCE, 
CE, MCE models [3] [4]. 

• This effect has been observed for the first 
time. 
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